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Abstract 

The nearly energy independent hadron emitter dimension r, measured in e^e~ annihi- 
lation in the energy range 10 to 91 GeV via the Bose-Einstein correlation of two identical 
charged pions, is shown to be well accounted for by choosing the hadron jets as inde- 
pendent pion sources. To this end the known normalised factorial cumulant moments 
dependence on particle sources is adapted to the Bose-Einstein correlation formalism to 
yield a relation between r and these sources. This approach is also able to account for the 
measured r values obtained for the decays into two and three hadron jets. Finally the 
estimated r value of the hadronic T(9.46) decay via three gluons is expected to be higher 
by about 6 to 11% over that predicted for its one photon hadronic decay mode. 
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1 Introduction 



The two-particle intensity interferometry has been proposed by Hanbury-Brown and Twiss [1] 
to correlate radio waves arriving from outer space for the purpose of measuring the angular 
diameter of astronomical objects. This method, which here is referred to as Bose-Einstein 
correlations (BEC), has been further extended by Goldhaber ct al., [2] to the system of identical 
hadron pairs produced in two-particle collisions. Taking the hadron emitter to be a sphere of 
radius r with a Gaussian distribution, the correlation function C(pi,p2) of pairs of identical 
bosons with momenta pi and p2 can be parametrised by 

C{Q) = 1 + Xe-'""^" , (1) 

where = —{pi — P2Y ■ The A parameter, usually referred to as the chaoticity parameter, 
is a measure of the incoherence level of the hadrons source. However, the estimated A values 
obtained from BEC analyses are also affected by the purity level of the identical hadron pairs 
sample taken from the experimental data. In recent years the BEC interferometry of pairs of 
charged pions is extensively used in order to investigate the underlying dynamical processes 
occurring in heavy ion collisions at ultra high energies [3]. In parallel, the BEC analysis was and 
is also applied to e^e~ annihilations at different centre-of-mass energies ^/s, and as a function 
of the outgoing hadron mass, in order to estimate the emitter dimension [4]. Lately a special 
attention has been given to the role, if any, of these BEC in the reaction e+e~ — > W^W~ — > 
hadrons in order to estimate their possible effect on the W mass determination [5] . 

An early compilation of the r^±7r± values measured in heavy ions collisions [6] is seen in 
Fig. 1 to be rather well described by the straight line representing the relation r = 1.2A^/^ fm, 
the known nucleus radius dependence on the atomic number A. This increase is attributed to 
the rise of the number of hadron sources, such as the nucleon-nucleon interaction, which occur 
as the atomic number of the colliding ions increases [3] . At the same time the BEC deduced r 
values, extracted from e'^e~ annihilations leading mainly to pions, seem to be rather indepen- 
dent of the centre of mass energy of the colliding electrons. This is illustrated in Fig. 2 where 
the small deviations from a constant r value are attributed to the variation in the experimental 
procedures adopted in each of the BEC analyses [4,7]. Additional information concerning the 
properties of the BEC extracted r dimension is coming from the hadronic decays where it 
is found that the r value increases with number of hadron jets and charged pion multiplicity [8]. 

Here it is shown that the observed energy independent behaviour of r in e+e" hadrons 
annihilations is consistent with the assumption that only hadrons emerging from the same final 
state hadron jet can be correlated. In addition, this approach leads also to an r increase of 
some 10% of the three-jet hadronic decays over that expected for two-jets decay events. 
Finally the hadron jet source approach is applied to the hadronic T(9.46) decays via the three 
gluons and one photon processes. 
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Figure 1: A compilation of r„±^± extracted from BEC analyses of identical charged pion-pairs pro- 
duced in heavy ion collisions shown as a function of A^/^ where A is the atomic number of the 
projectile. The values are taken from Ref. [6] and whenever more than one value was given for the 
same projectile the plotted value is their average. The line represents the relation r = 1.2A^/^ fm. 



2 The hadron source approach in cumulants and BEC 



A direct way to extract the genuine dynamical correlations of two and more hadrons, produced 
in high energy reactions, is offered by the so called bin-averaged normalised factorial cumulant 
moments Kg, here denoted simply as cumulants, where q is the size of the hadron group and 
is equal e.g. to 2 for a pion-pair system. This method, which was first proposed in [9] and 
described in details in [10], has been applied in the study of multi-particle dynamical fluctu- 
ations. In these investigations, also referred to as intermittency analyses, the dependence of 
the cumulants values on the average multiplicity has been proposed to be the consequence of a 
mixing of several well defined hadron emission sources [11-15]. 

In this emission source approach, correlations exist only between particles emitted from 
the same source so that the correlation strength is reduced whenever particles are grouped 
from more than one source. For a quantitative evaluation of this reduction we follow here 
closely Rcf. [15] in its treatment of the one dimensional (rapidity) cumulants K^'", which are 
averaged over M bins with the condition that M > 10. If one denotes the average value of a 
cumulant determined from one hadron source as K^"^^^ then the cumulant K^''^''\ calculated 
for all possible q-groups from a state k of more than one source, will be diluted by a factor D^'''> 
so that 

J^av{k) ^ ^(k) j^av{l) ^ (2) 
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Figure 2: Values of r^±7r± obtained from Bose-Einstein correlation analyses of hadrons produced in 
e'^e" annihilations [4,7]) plotted against ^/s. The shown band is the r^±7r± behaviour expected for 
a constant A value from the hadron source approach, normalised at 91 GeV to the r value obtained 
by the ALEPH collaboration. The continuous and dashed lines correspond respectively to the Xmax 
values of 0.99 and 0.95 used in calculating {Dqqg) via Eq. (12). 



where D^'^^ is equal to 



pG 



q ' ^ q 

Here denotes the number of q-particle groups, e.g. pairs of pions, emerging from the same 
source while P^'-^ stands for the number of all possible combinations of q-particlc groups which 
emerge from at least two sources. In Ref. [15] it has further been shown that in the case of S 
identical sources, each having the same charged multiplicity n, one has for g = 2 that 



This expression for = 2 can be generahsed for g > 1 to yield Dq 1/ S'^'^ . For the present 
work it is important to realise that even in the case where the number of sources remains the 
same, but the condition of equal multiplicity is lifted, the dilution factor Dq may change when 
the over-all multiplicity and/or its division between the sources is modified. 

In the case of two particle correlations, which obviously are genuine, one is able to relate [16] 
the normalised cumulant K2 and the BEC correlation function C(pi,p2), namely: 



C(Pl,P2) - 1 = i^2(Pl,P2) 
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(5) 



A comparison between this expression and Eq. (1) affords the possibihty to utihse the depen- 
dence of the averaged normahsed cumulant on the hadron sources and their charged multiphcity 
configuration to the BEC correlation term Xe^^ ^ average over the Q variable. Specifically, let 
us consider the ratio of the normalised average cumulants K^^^^ / K2^^^\ where the superscripts 
m and k denote two different hadron sources configurations. Thus one has 

^0 ' Jo XkVm 

Prom Eq. (2) follows that = D^^ / o'i^ so that one finally obtains a simple 

relation between the hadron emitter radii and the dilution factors, namely 



Xm D2 



(k) 

Tk . (6) 



If wc set the final state k to be that emerging from a single source and assign to it the 
value k = 1, then D^^ serves as a yardstick against which dilution factors of other final states 
are measured. Thus for identical sources having the same A and average charged multiplicity, 
Eq. (6) reduces with the help of relation (4) to = Sri where ri is the emitter radius 
of a single hadron source. Hence r increases linearly with S, the number of hadron sources. 
Although it seems a priori naturally to choose for the heavy ions reactions the nucleon-nucleon 
collision as the basic single hadron source, the actual description of nuclei reaction is much 
more complicated mainly due to the secondary interactions within the nucleus. 



3 r7r±7r± versus ^/s in e+e annihilation 

In the following we propose to study the dependence of r^±^± on y^, as measured in the 
e~^e~ hadrons annihilation, in terms of the hadron sources given by the quark and gluon 
jets. Experimentally, for the grouping of the outgoing hadrons into jets there exist several 
methods which differ somewhat in their results. Among the leading ones is the known Durham 
algorithm [17], which has been applied to the hadronic decay. The outcome of such a hadron 
jet analysis, as obtained by the OPAL collaboration [18], is shown in Pig. 3 as a function of 
the jet resolution parameter ijcut- As can be seen, already for values of ijcut ^ 0.01, the data is 
predominantly limited by three hadron jets. In the following we will assume that at s/s < 91 
GeV the e'^e~ — > hadrons annihilation is governed by not more than the three-jet configuration. 

In the first next-to-leading process e"^e~ — > qqg, the differential cross section is given by 
(see e.g. Ref. [19]): 

adxidx2 ~ Stt (l-a;i)(l-a;2) ~ Stt '^^^^^^i' ^^j ■ 10 

Here Xi = lE^j ^ and X2 = 2Eg/y/s. where Eg and Eg are respectively the energies of 
the outgoing quark and anti-quark jets. At this point it is important to realise that the term 
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Figure 3: The relative measured rate R of the produced hadron jets in the decay, taken from 
reference [18], as a function of the jet resolution parameter ycut using the Durham algorithm. The 
data which is corrected for all the OPAL detector deficiencies, but not for hadronisation effects, is 
compared with two Monte Carlo models. 

Wqqg^xi, X2) is equivalent to the probability of finding a gluon with the energy Eg = y/s—Eg—Eg 
which here is utilised for the calculation of the dilution factors. The integration region asso- 
ciated to this differential cross section is defined by < Xi, 0:2 < 1 and Xi + X2 > 1. Further 
to note is that the expression given by Eq. (7) diverges when either the gluon is coUinear with 
one of the outgoing quarks or the gluon momentum approaches zero. Some methods to handle 
in the integration these singularities are discussed e.g. in references [19,20]. Here we deal with 
this singularity by introducing to the Xi and X2 variables an upper limit of x^ax — 0.99. To 
estimate the sensitivity of our calculations to this upper limit cutoff choice we also represent 
our analysis results using Xmax — 0.95. 

For the average charged multiplicity of the quark and anti-quark jets, in the process e'^e~ — > 



qqg, it was found that its parametrisation is best given in terms of an appropriate scale defined 
as [21] 



where Eq (Eq) is the energy of the quark (anti-quark) jet and 6qg {9qg) is the opening angle 
between the quark (anti-quark) jet and the gluon jet. The effective scale for the gluon is then 



For the dependence of the multiplicity on the scale Qjet we use the parametrisation given 




and 






(9) 
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in Refs. [19,21], namely: 

{NqiQjet)) = ao + ailnQjet + a2{lnQjety 



(10) 



and 

(NgiQjet)) = Ro + Rl{N,{Qj,t)) : (11) 

where Oq = 2.74 ± 0.07, Oi = 1.71 ± 0.05, 02 = -0.05 ± 0.07, Rq = -7.27 ± 0.52 and 
Ri — 2.27 ± 0.07. With the help of these two parametrisations one can calculate the 2-pion 
dilution factor Dqqg{xi,X2), as defined in Eq. (3), for every set of xi, X2 and Xg = 2 — xi — X2 
values assuming for simplicity that each jet charged multiplicity is equally divided into positive 
and negative pions. Thus the average 2-pion dilution factor (Dggg) , at a given e'^e" centre-of- 
mass energy, is equal to 

r^max f'Xmax I P^max f'^max 

{Dqqg) = j dX2 Dqqg{xi,X2)Wqqg{Xi,X2)dXi/ / dX2 Wqqg{Xi,X2)dXi . 

Jo JI—X2 JO JI—X2 

(12) 

The average dilution factors {Dqqg) were thus calculated at several ^/s values in the range of 
10 to 91 GeV using the two Xmax values of 0.99 and 0.95. Our results for the expected r^±7r± 
dependence on energy, using Eq. (6) and assuming a constant A value, are represented in Fig. 2 
by a band chosen to be normalised to the r value measured by ALEPH [22] at 91 GeV so as 
to cover also r measurements at lower energies. The lower and upper lines of the band, which 
correspond respectively to the x„,,a:r vahies of 0.99 and 0.95, are seen to be very similar in their 
shape. As can be seen, the expected dependence of r on the energy is weak and consistent with 
the experimental observations in the energy range from 10 GeV up to the r values measured at 
the Z'^ mass energy by the ALEPH and DELPHI collaborations. At the same time, in the frame 
work of the hadron source approach, no relation can be achieved between the reported r values 
of L3 and OPAL at = 91 GeV and those measured at lower energies. Finally important to 
note is that whereas in the description of the jet multiplicity via Eqs. (10) and (11) the use of 
Qjet scale is mandatory, [21], the average calculated dilution factor (Dqqg) differs, in the 10 to 
91 GeV energy range, only by about 2% when the Qjet scale is replaced by the jet energy Ejet- 

3.1 The r7r±7r± value of the decay into two and three jets 

The dependence of the emitter size on the number of hadron jets in the Z° decays has been 
studied by the OPAL collaboration [8] which found that the r value of the 3-jet events has a 
value which is higher by about 10% than that obtained for the 2-jet events. From the same 
study OPAL also reported that r increases with the observed charged multiplicity while the A 
parameter decreases. As mentioned before, the experimental separation between 2-jet and 3-jet 
events depends on the method chosen and its parameters values and therefore its reproduction 
with Eq. (7) is inaccessible. On the other hand, transition from a 2-jet to a 3-jet decay con- 
figuration can be characterised, in the framework of the hadron source model, by an increase 
of A''^, the number of charged hadrons associated with the gluon jet. The expected ratio of 
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Figure 4: The expected ratio r{Ng)/r{Z^) of the hadronic decay as a function of Ng, the gluon 
jet charged multiphcity, as calculated from the hadron source approach. 

r{Ng)/r{Z^), as obtained from the model using x^ax — 0.99, is shown in Fig. 4 as a function of 
Ng, where r(Z°) is the dimension value obtained from the BEC analysis of the total Z° final 
state hadrons. As seen, this ratio increases by about 10% from Ng of about three to values above 
seven where the ratio settles at a level of about ~1.02. Similar behaviour of r{Ng)/r{Z^) is also 
seen when the value Xmax — 0.95 is used in the calculations. Since r(Z°) lies between the r val- 
ues of the 2 and 3-jet events, the model expectation is consistent with the experimental findings. 

For the comparison of the OPAL results for the r dependence on multiphcity at the Z^ decay 
[8], with the hadronic jet source approach one requires above all a parametrisation of the hadron 
multiplicity spectrum as a function of Qjet in addition to its average value {N{Qjet)). For the 
gluon jet this requirement can be achieved by using its hadron multiplicity distributions given 
in [23] for several Qjet values. As for the quark jet, one can obtain the needed parametrisation 
by using the known average multiplicity versus Qjet and by assuming that the shape of the 
multiplicity distribution is independent of the jet energy. With these parametrisations we find 
that for a constant A the dimension r indeed increases with the true multiplicity. Nevertheless, 
it is not surprising that we are unable to reproduce the experimental findings in all their details 
and that for the following reasons. Firstly, the measured multiplicity given in [8] is the observed 
one, which depends on the experimental setup and analysis procedures, whereas the hadron 
jets parametrisations refer to the true multiplicity. Secondly, the measured chaoticity A is 
found to decrease rather steeply with the multiplicity which may indicate an increase in the 
experimental data contamination as the number of outgoing particles increases. Thirdly, the 
present study has limited itself to the case where the decay configurations of more than three 
hadron jets can be neglected. While this assumption is supported by the OPAL analysis [18] 
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even at ^/s — 91 GeV, as long as the BEG analysis embraces the whole final state hadrons, 
it may not be anymore valid when only the higher multiplicity part of the hadronic Z° decay 
events is selected for the correlation study. 



3.2 Quarkonium decay into hadrons 

Quarkonia, like the J/iIj{3.1) and the T(9.46), decay into hadrons predominantly via the three 
gluon {3g) mode, also referred to as the direct decay mode, which leads to three hadron jets. In 
addition, they also decay via the one photon, or vacuum polarisation, mode which is identical 
in its hadron jet configuration to that of the e+e" —>■ j/Z* qqg background underneath. 
The J/^(3.1) mass is too low to be analysed here in terms of its hadron sources, mainly due 
to the uncertainty in the gluon and quark jet parametrisations below ~ 1 — 2 GeV. As for 
the T(9.46), its mass may be considered just heavy enough to render, according to Eq. (6), a 
reliable estimate of the ratio 

/n \ \„ 

(13) 



i-^Sg) \qg 

To calculate the dilution factor of the three gluon-jet configuration we consider the differential 
T(9.46) direct decay rate given by [24] 

TT-:r^ = -^-^ W{x^,X2,Xs) , (14) 
r^g dXidX2 TT^ — 9 

where Xj = 2Eg./mr- Here W{xi,X2,X3) is given by [24] 

W{xi, X2, xs) = [xl{l - xif + xl{l - X2f + xl{l - x^f\ (15) 

X'^X2X^ 

with the condition that xi -\- X2 -\- x^i — 2. Whereas for the calculation of the average dilution 
factor {Dqqg) one has a well defined Qj^t scale for the quarks and gluon jets as formulated in 
Eqs. (10) and (11), not so for the three gluon configuration. However, as mentioned before, the 
average dilution factors do change only by ~ 2% when the Qjet scale is replaced by Ej^-t. This 
being the case we proceeded to calculated r^g/r-qqg using the jet energy Ej^t rather than the Qj^t 
scale while incorporating a 2% uncertainty into our final result. The evaluation of the average 
dilution factor (-Dsp) is then carried out via Eq. (12) where Dqqg{xi,X2) and Wqqg{xi,X2) are 
replaced respectively by D^g{xi,X2) and Wsg{xi,X2). From these calculations and those carried 
out to obtain (Dqqg) at ^/s = 9.46 GeV, we find that the ratio (Dqqg) / (D-^g) is lying in the 
range of 1.06 to 1.11 corresponding to the integration cutoff values Xmax = 0.99 and 0.95. 

An early BEG analysis of the GLEG collaboration [25] yielded the values of r(T — ^ 3g) = 
0.99±0.14 fm and r(T qqg) = 0.86±0.15 fm so that their ratio is equal to 1.15±0.26. A more 
recent BEG measurements of e~^e~ — > hadrons at the T mass energy, described in Ref. [26], 
quotes the values r(T 3g) = 0.73 ± 0.10 ± 0.04 fm and r(T qqg) = 0.83 ± 0.22 ± 0.05 
fm, from which follows that r(T — > 3g)/r{T qqg) = 0.88 ± 0.28. To compare these results 
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with the expectation of the hadron source approach one needs to know the true X^g and Xqqg 
values. However as mentioned before, these are difficult to estimate due to the fact that 
they are sensitive to the purity level of the experimental data sample. Inasmuch that X^g ~ 
Xqgg, an assumption supported by the analysis results reported in [26], the two independent 
measurements of r^g/rggg are consistent within one standard deviation with the expectation of 
the hadron source approach. 

4 Summary 

The dependence of the Bose-Einstein correlation of two identical bosons, produced in e+e~ 
annihilations, on the configuration of the hadron sources has been taken over from the known 
cumulants dependence on these sources. Assuming that correlation exists only for hadrons 
emerging from the same source, leads to a BEC dilution whenever the hadron-pair emerges 
from two different sources. A relation between the amount of dilution and the BEC extracted 
dimension r is formulated. Specifically, we identify in the e^e~ annihilation the outgoing parti- 
cle jets as the hadron sources which up to centre-of-mass energy of ^Js — 91 GeV are essentially 
limited to three. 

The comparison of the hadron source approach to the experimental findings is somewhat 
hindered by the need to know the genuine values of the chaoticity A parameters, and in some of 
the cases also the distribution of the true multiplicity, both as a function of the jet energy. In 
addition, the extension of the hadron source approach to e^e^ annihilations at energies above 
91 GeV seems to require the inclusion of configurations higher than three hadron jet sources 
which are not dealt with here. 

Insofar that A remains the same over the whole energy range from 10 to 91 GeV, the hadron 
source approach expects the BEC extracted r value to be only slightly dependent on energy 
in accordance with the experimental observations. As for the hadronic decay, the seen rise 
of r by some 10% when shifting from the two to the three hadron-jet configuration, is also 
expected from the hadron source approach. Finally the model is consistent within errors with 
the experimental findings for the ratio between the r values obtained for the three gluons and 
to the one photon hadronic decay modes of the T(9.46). 
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